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SUMMARY We give a centralized deterministic algorithm 
for constructing linear network error-correcting codes that attain 
the Singleton bound of network error-correcting codes. The pro- 
posed algorithm is based on the algorithm by Jaggi et al. We 
give estimates on the time complexity and the required symbol 
size of the proposed algorithm. We also estimate the probability 
of a random choice of local encoding vectors by all intermedi- 
ate nodes giving a network error-correcting codes attaining the 
Singleton bound. We also clarify the relationship between the ro- 
bust network coding and the network error-correcting codes with 
known locations of errors. 

key words: error correction, MDS code, network coding, ran- 
dom network coding, Singleton bound 

1. Introduction 

Ahlswede et al. [1] proposed the notion of network cod- 
ing that multicasts data from a single sender to multi- 
ple receivers at a rate at which the ordinary store and 
forward routing cannot multicast the data. Such high 
rate multicast becomes feasible by allowing intermedi- 
ate nodes to encode and decode the data. A sender is 
usually called a source and a receiver is called a sink. A 
network coding is said to be linear if every intermediate 
node outputs a linear combination of its inputs [10] . 

A study of network coding usually assumes that 
an error does not occur in networks. Recently, Cai and 
Yeung [2], [13] considered errors in network coding, and 
proposed the network error correcting codes that allow 
sinks to recover the information even when errors oc- 
cur on intermediate edges in the network. After formu- 
lating the network error correction, they proposed the 
lower and upper bounds on the number of messages 
in a network a-error correcting code, and one of their 
upper bound was a natural generalization of the Sin- 
gleton bound for the ordinary error-correcting codes. 
Recently, Zhang [14] and Yang et al. [11] independently 
observed that the Singleton bound can be refined. We 
note that the problem formulation in [2], [13] was later 
independently presented in [4] . (The proceedings paper 
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of [2], [13] appeared in 2002.) 

Cai and Yeung mostly considered the case that in- 
termediate nodes perform only simple encoding and de- 
coding without delay, such as computing the output of 
the node as a linear combination of its inputs, and the 
sinks perform complex decoding computation. The net- 
work error correcting codes can avoid introducing de- 
coding computation and delay into intermediate nodes, 
which is the advantage over use of ordinary error cor- 
recting codes between nodes. 

Note that a similar type of network failure in a 
slightly different context was considered in [7, Sect. V] 
and [6, Sect. VI] in which every sink is assumed to know 
the set of failed edges and failed edges are assumed 
to emit zero symbols. Network error correction does 
not assume the knowledge of edges causing errors, and 
the problem formulation is different from [6] , [7] . Note 
also that Kurihara [8] considered the different notion of 
robustness. In his paper, he considered network coding 
that allows sinks to recover partial information with 
edge failures. 

For the construction of the network error- 
correcting codes, Jaggi et al. [5] proposed a random- 
ized construction that uses coding among different time 
intervals. Their method produces codes attains the 
Singleton bound with high probability with sufficiently 
long block length, where the block length refers to the 
number of time intervals among which coding is done. 
It is desirable to have a network error-correcting code 
that does not code among different time intervals and 
thus does not introduce delay. Concurrently to this pa- 
per, Yang et al. [11] proposed an explicit construction 
algorithm that produces codes attaining the refined Sin- 
gleton bound. The idea in [11] is similar to this paper 
in the sense that they also regard errors as information 
from the source and add extra components in the global 
encoding vectors corresponding to errors. 

In this paper, we give a deterministic and cen- 
tralized algorithm that constructs a network error- 
correcting code that attains the Singleton bound of 
network error-correcting codes obtained in [13]. We 
also give a relationship between the success probability 
and the field size for successful construction of network 
error-correcting codes when intermediate nodes choose 
their encoding coefficients randomly and independently 
The proposed algorithms are based on [6] . Our network 
error-correcting codes make multicast robust to errors 
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without introducing delay in the transmission, which 
is very attractive to delay sensitive multicast applica- 
tions, such as multicast of video or audio. Our method 
is also useful for cryptographic applications, because it 
can tolerate modification and deletion of data by an 
adversary. 

This paper is organized as follows. Section 2 in- 
troduces notations and the model of errors. Section 3 
proposes an algorithm for constructing network error- 
correcting codes attaining the Singleton bound. Section 
4 shows how to modify the algorithm in Sect. 3 to attain 
the refined Singleton bound, the success probability 
of the random construction of network error-correcting 
codes, and the relationship between the robust network 
coding [6] , [7] and the network error-correcting codes 
with known locations of errors [12]. Section 5 gives 
concluding remarks. 

2. Preliminary 

2.1 Basic notations 

We consider an acyclic directed graph G = (V, E) with 
possible parallel edges of unit capacity. V 3 s denotes 
the source and V D T denotes the set of sinks. Let n 
be the smallest min-cut separating s from any t £ T 
throughout this paper. For v G V, T + (v) (resp. r~(v)) 
denotes the set of edges leaving (resp. reaching) the 
node v, and start(e) (resp. end(e)) denotes the node at 
which the edge e starts (resp. ends). 

We consider linear coding over a finite field F q with 
q elements. The source s gets k (< n) input symbols 
from F q . The symbol y(e) G F q carried by an edge e 
is a linear combination of the symbols carried by the 
edges entering start(e). The local encoding vector m e : 
r _ (start (e)) — > F q determines the coefficients of this 
linear combination, that is, 

y(e) = m e (e')y(e). 

e'er-(start(e)) 

In this paper, a nonsink node performs only the 
computation of linear combination of its inputs, and 
they do not correct errors. An error is assumed to occur 
always at an edge. When an error occurs at an edge e, 
the symbol received by end(e) is different from one sent 
by start(e), and end(e) computes its outputs as if there 
was no error at e. The error value at an edge e is defined 
by the received symbol minus the transmitted symbol 
at e. Note that we express a failure of a node ueKin 
a real network as errors on edges in T + (v) in our model. 
The number of errors is the number of edges at which 
errors occur. A network code is said to correct a errors 
if every sink can recover the original information sent 
by the source when a or less errors occur at arbitrary 
edges. We call the recovery of information by a sink 
decoding. 

We represent errors occurred in the whole network 



by a vector e in Fj, , where \E\ denotes the number 
of elements in E. Fix some total ordering in E, and 
enumeration of the error values gives e. 

Regarding on the number of messages in a network 
en-error correcting code, Cai and Yeung obtained the 
following result. 

Proposition 1: [13] The number M of messages in a 
network a-error correcting code, not necessarily linear, 
is upper bounded by 

M < q n - 2a . 

Very recently, Zhang [14] and Yang et al. [11] ob- 
served that the above proposition can be refined as fol- 
lows. 

Proposition 2: [11], [14] Let n t be the min-cut from 
the source s to a sink t. If the sink t can correct any 
a t errors then the number M of messages in the net- 
work correcting code, not necessarily linear, is upper 
bounded by 

M < q n t- 2a \ 

2.2 Jaggi et al.'s algorithm for construction of an or- 
dinary network code 

In this subsection, we review Jaggi et al.'s algorithm [6] 
for construction of an ordinary network coding. The 
proposed algorithm uses a modified version of their al- 
gorithm. 

Since linear coding is used, the information carried 
by an edge e is a linear combination of k information 
symbols in F q . We can characterize the effect of all the 
local encoding vectors on an edge e independently of 
a concrete k information symbols using global encoding 
vectors b(e) € F*. When the information from the 

source is i G F q , the transmitted symbol on an edge 

e is equal to the inner product of i and 6(e). In order 
to decide the encoding at the source node s, we have 
to introduce an imaginary source s' and k edges of unit 
capacity from s' to s. We regard that s' sends k symbols 
to s over k edges. 

We initially computes an s'-t flow /* of magnitude 
k for each t G T and decomposes this flow into k edge 
disjoint paths from s' to t. If an edge e is on some flow 
path W from s' to t, let /^_(e) denote the predeces- 
sor edge of the edge e on the path W. Jaggi et al.'s 
algorithm steps through the nodes v G V in a topo- 
logical order induced by the directed graph G. This 
ensures that the global encoding vectors of all edges 
reaching v are known when the local encoding vectors 
of the edges leaving v are determined. The algorithm 
defines the coefficients of m e for one edge e G T + (v) 
after the other. There might be multiple flow paths 
to different sinks through an edge e. Let T(e) de- 
note the set of sinks using e in some flow /* and let 
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P(e) = {/f_(e) | t £ T(e)} denote the set of predeces- 
sors edges of e in some flow path. The value is chosen 
for m e (e') with edges e' £ P(e). 

We introduce two algorithmic variables B t and Ct 
that are updated by Jaggi et al.'s algorithm. C t con- 
tains one edge from each path in /*, namely the edge 
whose global encoding vector was defined most recently 
in the path. B t = {b(e) \ e £ Ct} is updated when Ct 
is updated. The algorithm determines m e so that for 
all t £ T, B t is linearly independent. 

After finishing the algorithm, every sink can de- 
code the original information because B t is linearly in- 
dependent. 

3. Construction algorithm 

We shall propose an algorithm constructing a network 
a-error correcting code carrying k information symbols 
in F q with n — k > 2a, which is equivalent to the Single- 
ton bound (Proposition 1). The proposed construction 
is based on [6] . We assume that the size of alphabet F q 
satisfies 



q>\T\- 




Definition 3: For the original information i £ F^ 

and the error e £ f[ E \ let (j>t{i,e) £ Fq F ^ be the 
vector of symbols carried by the input edges to t. 

Lemma 4: If a sink t can decode the original infor- 
mation i with any 2a or less errors whose locations are 
known to the sink t, then the sink t can decode the orig- 
inal information with any a or less errors without the 
knowledge of the error locations under the assumption 
that the number of errors is < a. 

Note that errors with known locations are called era- 
sures in [12] and the properties of erasures are also stud- 
ied in [12]. 

Proof. Denote the Hamming weight of a vector x by 
w(x). The assumption of the lemma implies that for 
any i ^ j and e with w(e) < 2a we have 

<M?,e)^<MJ,0). (2) 

Equation (2) implies that for any i ^ j and e\, ei with 
w(ei) < a and w{e2) < a we have 

<t>t{i,ei) ^ <M/,e 2 ), 

which guarantees that t can decode the original infor- 
mation under the assumption that the number of errors 
is < a by exhaustive search. □ 

Remark 5: The above lemma does not guarantee the 
existence of an efficient decoding algorithm. 

Fix F C E with \F\ = 2a. We shall show how to 




Fig. 1 Example of a network with imaginary nodes and edges. 
Nodes A and B are the imaginary nodes added in Step 2 and the 
dashed lines from s' to A and B represent the imaginary edges 
added in Step 2. See Example 6 for explanation. 

construct a network error-correcting code that allows 
every sink to decode the original information when the 
errors can occur only at F. We call F the error pat- 
tern. The following description is a condensed version 
of the proposed algorithm, which is equivalent to the 
full description with T = {F} in Fig. 2 on p. 5. 

1. Add the imaginary source s' and draw k edges from 
s' to s. 

2. Add an imaginary node v at the midpoint of each 
e £ F and add an edge of unit capacity from s' to 
each v. 

3. For each sink t, do the following: 

a. Draw as many edge disjoint paths from s' to 
t passing through the imaginary edges added 
at Step 2 as possible. Let mf (< 2a) be the 
number of paths. 

b. Draw k edge disjoint paths passing through s 
that are also edge disjoint from the mf paths 
drawn in the previous step. 

4. Execute the algorithm by Jaggi et al. with 
^ teT (fc + mf ) edge disjoint paths constructed in 
Step 3. 

Example 6: In Fig. 1, we give an example of addition 
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of imaginary nodes and edges. The network structure 
in Fig. 1 is taken from [3, Fig. 2]. Nodes A and B are 
the imaginary nodes added in Step 2 and the dashed 
lines from s' to A and B represent the imaginary edges 
added in Step 2. 

The min-cut from s to every sink is 4 in the original 
network. The set F of edges with errors consists of the 
edge from s to node 1 and the edge from node 1 to node 
5. 

We denote a path by enumerating nodes on the 
path. In Step 3a for t\ we can find two edge dis- 
joint paths, namely (s' , A, 1, t\) and (s' , B, 5, 8, t\). On 
the other hand, in Step 3a for t 2 , we can find only 
one edge disjoint path, namely (s', A, 1, B, 5, 8, t 2 ) or 
(s', B, 5, 8, t 2 ). Therefore m£ = 2 while m£ = 1. 

In Step 3b for t\, we find two edge disjoint paths 
as (s', s, 3, 6, 9, t x ) and (s', s, 4, 7, 10, t x ). In Step 3b for 
t 2 , we find three edge disjoint paths as (s', s, 2, 6, 9, £2), 
(s', s, 3, 7, 10, i 2 )> an d (s' , 4, i 2 )- We can use arbitrary 
two paths among the three paths. In either case, we can 
find n — mf paths in Step 3b. □ 

In Step 3b, we can guarantee the existence of k 
paths as follows: Suppose that edges in mf paths 
used in Step 3a are removed from the original network 
(V, E) . Then the min-cut from s to a sink t in the orig- 
inal network (V, E) is at least n — mf, which is larger 
than or equal to k. 

In Step 4 we use the algorithm by Jaggi et al. as 
if the imaginary source s' sent information on the a 
imaginary edges added in Step 2. We denote by Bf 
the set B t of global encoding vectors for k + mf edge 
disjoint paths. Bf consists of k + mf vectors of length 
k + 2a. We require that every sink t is able to decode k 
information symbols, while t may be unable to decode 
2a error symbols in general because mf < 2a. 

There are always two edges end at the added imag- 
inary node v and one edge starts from v in Step 2. Since 
v is imaginary, we cannot choose local encoding vectors 
at v. Therefore, in Step 4, all components in the local 
encoding vector at v must be selected to 1, which keeps 
B t linearly independent. The reason is as follows: Let 
e be the edge from s' to v added in Step 2. The global 
encoding vector of e is of the form 

(o j -\i,o n - j ), 

that is, it has only 1 at the j-th component. All other 
global encoding vectors in Bf have zero at the j-th 
component, since they are not in downstream of e when 
we choose local encoding vectors at v. Therefore, the 
added imaginary node v does not interfere with the 
execution of Jaggi et al.'s algorithm. 

Observe also that q > \T\ guarantees the successful 
execution of the algorithm as with the original version 
of Jaggi et al.'s algorithm. 

We shall show how each sink t can decode the origi- 
nal information sent from the source s. After executing 



Step 4 we have decided all the local encoding vectors in 
the original network (V,E). Consider the three linear 
spaces defined by 

V x = {&(i,g) |?eF£,eeFfl}, 

V 2 = {MhO)\^F k q }, 

V 3 = {<M0,e) |e£Fl £ l}, 

where components in e corresponding to E\F are zero, 
and <fit is as defined in Definition 3. We consider V\, 
V2, and V3 in the original network (V, E) without added 
imaginary nodes and edges. Then we have 

^1 = ^2 + ^3, dimV 2 <fc. (3) 

Since we keep Bf linearly independent, 

dim^i > k + mf. (4) 

Since the maximum number of edge disjoint paths pass- 
ing through the imaginary edges added in Step 2 is mf, 
we have 

dim V 3 <mf. (5) 
Equations (3-5) imply 

dim Vi = k + mf , 

dimV 2 = k, (6) 
dimV3 = mf, 
dimV 2 r\V 3 = 0. (7) 

The number of nonzero components in (fi t (i,e) is k+mf 
and the number of unknowns in <fit(i, e) is k + 2a, which 
can be larger than k+mf . However, by Eq. (7), the sink 
t can compute (fi t (i,0) from (fit(i,e) as follows: Write 
(fit{i,e) as u + v such that u G V 2 and v E V 3 . By Eq. (7) 
u and v are uniquely determined [9, p. 19, Theorem 4.1]. 
We have u = 4>t(i,0) and the effect of errors is removed. 
The sink t can also compute the original information i 
from (fi t (i,0) by Eq. (6). 

We shall describe how to construct a network error- 
correcting code that can correct errors in any edge set 
F C E with |F| = 2a. Let T = {F C E : \F\ = 2a}. 
The idea in this paragraph is almost the same as the 
construction of the robust network coding in [6, Sect. 
VI] . Recall that Bf is the set of global encoding vectors 
on edge disjoint paths to a sink t with an edge set F 
of errors. Execute Jaggi et al.'s algorithm keeping Bf 
linearly independent for all t € T and all F G T ' . Then 
every sink t can decode the original information with 
the knowledge of the edge set F on which errors actually 
occur. As in [6, Sect. VI], 

,>m-i^i = iTi(fJ) 

guarantees the successful execution of the algorithm. 
We present a pseudo programming code of the 
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(* Initialization *) 

Added imaginary node s' and edges ei, . . . , efc from s' to s. 0(k) 
foreach error pattern F £ T do 
Initialize global encoding vector 

b F (e x ) = (O*- 1 , 1, fe + 2 «- i ) £ F£+ 2< \ 0((fc + 2a) 2 ) 

foreach edge e 6 F do 

Add an imaginary node v at the midpoint of e 6 F.O(l) 
Divide e into an edge to v and an edge from v. O(l) 
Draw an imaginary edge from s' to u. (*) O(l) 

endforeach 

foreach sink i£Tdo 

Draw as many edge disjoint paths from s' to t as possible 
passing through the edge added in (*). 

0(2a(\E\ + k + 4a)) 
Draw k edge disjoint path from s' to t passing through 
s and also disjoint from paths made in the previous step. 

0{k(\E\ + k + 4a)) 
Initialize the basis Bf = {b F (ei) | ej is on a path to t}. 

0((k + 2a) 2 ) 

endforeach 
endforeach 

(* Main loop *) 

foreach edge e g {Jp e jr Ef \ { e li ■ • ■ > e fc} m a topological order do 

if start (e) e V then 

Choose a linear combination b F '(e) = ^2 P eP F (e) me (p)^(p) 
such that B F remains linearly independent for all t 
and F by the method in [6, Sect. III.B]. (**) 

else 

m e (p) = 1 for all p e P F (e) and b F (e) = £ pepF(e) ?(p). 

0(fc + 2a) 

endif 
endforeach 

return {m e (-) | start(e) S V}. 

Fig. 2 Construction algorithm for a network o-error correcting 
code. The rightmost O(-) indicates the time complexity executing 
the step. 



proposed algorithm in Fig. 2. In order to present 
a detailed description, we introduce new notations. 
Gp = (Vp,E F ) denotes the network with added imag- 
inary nodes and edges in Steps 1 and 2 with the error 
pattern F C E. Let / t,F be the flow established in 
Steps 3a and 3b in G F - Let /^i F (e) denote the set of 
predecessor edges of the edge e in a flow path in f tF . 
Let T F {e) denote the set of sinks using e in some flow 
f- F and let P F {e) = {/1(e) | t G T(e)}. 

We shall analyze the time complexity of the pro- 
posed algorithm in Fig. 2. As in [6] we assume that 
any arithmetic in the finite field is O(l) regardless of 
the field size. First we analyze that of the initializa- 
tion part. Observe that \E F \ = \E\ + k + 2\F\ = 
\E\ + k + Aa because each edge in F adds two edges 
to E and there are k edges from s' to s. The most time 
consuming part in the initialization is construction of 
edge disjoint paths, whose overall time complexity is 
0((\E\ +k + Aa)\T\\T\{k + 2a)). 

Next we analyze the time complexity of the main 
loop. By [6, Proof of Lemma 8], the time complexity of 
choosing the local encoding vector m e (p) in Step (**) is 
0({\T\\T\) 2 {k+2a)) , which is the most time consuming 



part in the main loop. Choice of m e (p) is executed for 
\E\ edges starting from a real node in V. Thus, the 
time complexity of the main loop is Od-Eld^HTl) 2 ^-!- 
2a)), and the overall time complexity is 0(|.F||T|(fc + 
2a)[\E\ + k + Aa+\J : \\T\]). Note that \F\ = ( |F| ). 

A sink decodes the information by exhaustive 
search. Specifically the sink enumerates all the pos- 
sible information and all the possible errors for all 
F e T, then compares the resulting symbols on in- 
coming edges with the actual received symbols by the 
sink. The computation of the resulting symbols can be 
done by a matrix multiplication in 0((k + a) 2 ) time 
complexity. The number of possible information is q k 
and the number of possible errors is X)J=o (' ^'X^ - 
Thus, the time complexity of decoding by a sink is 

o(9 fe E; = o('?)(<7-i) j (fc+«) 2 )- 

4. Variants of the proposed method and its re- 
lation to the robust network coding 

We shall introduce two variants of the proposed method 
in this section. 

4.1 Attaining the refined Singleton bound 

Network error-correcting codes constructed by the pro- 
posed method attains the Singleton bound (Proposition 
1), while they do not necessarily attains the refined 
Singleton bound (Proposition 2). Yang ct al. [11] con- 
currently proposed a construction algorithm that pro- 
duces a code attaining the refined Singleton bound. In 
this subsection we modify the proposed method so that 
it can produce a code attaining the refined Singleton 
bound. 

Let n t be the min-cut from s to t, and suppose that 
the source s emits k symbols within unit time interval. 
A sink t can correct a errors if 2a < n t — k. Let J- t — 
{F C E : \F\=n t - k} and T = UtGT-^t- For fixed 
F e T and t G T, we cannot garuantee that there exists 
k edge disjoint paths in Step 3b. For such F, the sink 
t cannot decode information with errors occered at F. 
We exclude Bf with such (i, F) from the algorithm. 
Note that if < n t — k then there always exist k edge 
disjoint paths in Step 3b. 

In order to attain the refined Singleton bound we 
keep the linear independence of all bases in {Bf \ t G T, 
F G \F\ < n t - k} in Step (**) in Fig. 2. By the 
exactly same argument, we see that the produced code 
attains the refined Singleton bound. 

By almost the same argument as Sect. 3, we see 
that the modified proposed algorithm runs in time com- 
plcxity 0(\T\\T\{k + 2a max ) [\E\ + k + 4a max + |JT| |T|]), 
where a max = |_(max tg T n t — k) /2J . The required field 
size for successful execution of the algorithm is \T\ ■ 
and in this case \T\ depends on the structure of the 
network (V, E) . 
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Table 1 Comparison among the proposed methods and [5], [11]. Wc assumed that the 
min-cut is n for all t € T and k = n — 2a. I denotes the maximum of in-degrees of nodes. 





delay 


required field size for the 
success probability of code 
construction to be > 1 — 8 


time complexity of code construc- 
tion 


time complexity of decod- 
ing by sinks 


Figure 2 


none 




0(([^)\T\(k + 2a)[\E\ + k + A a + 


0(5*E°'-o( l f l )(9- 1 ) , '( fc + 
«) 2 ) 


Sect. 4.2 


none 


i^im (£')/« 


O(I) 


0(5*E°'-o( l f l )(9- 1 ) , '( fc + 
«) 2 ) 


Paper [11] 


none 


| T |(n +£ |- 2) 


O(\E\\T\ q k^ a =0 (tf)( q -iy) 


0(5 fc E;-o( l f')(9- 1 ) , '( fc + 


Paper [5] 


large 


not estimated 


O(I) 


0((n X delay)" 3 ) 



On the other hand, the time complexity of con- 
structing local encoding vectors by the method of Yang 
et al. [11] is 




and the required field size is 
v fn t + \E\ - 2\ 

The time complexity of the proposed algorithm can be 
smaller or larger depending on the network structure 
and q than Yang et al. [11]. The required field size of 
the proposed algorithm can also be smaller or larger 
depending on the network structure. However, for the 
special case n t = n for all t £ T, the required field size 
of the proposed method is smaller than Yang et al. [11]. 

4.2 Completely randomized construction 

By using the idea in the previous section, we can esti- 
mate the success probability of constructing a network 
error-correcting code by randomly choosing local en- 
coding vectors as follows. The idea behind its proof 
is almost the same as [6, Theorem 12]. Observe that 
the random choice of local encoding vectors completely 
remove the time complexity of selecting encoding vec- 
tors in the centralized manner at the expense of larger 
required field size q. 

Proposition 7: Suppose that the source s transmits 
k symbols within unit time interval, and let T — {F C 
E : \F\ — 2a} be the set of edges on which errors 
can occur. Suppose also that local encoding vector co- 
efficients are generated at random independently and 
uniformly over F q . With this network error-correcting 
code, all sinks can correct errors in any edge set F E T 
with probability at least 1 - S if q > \E\\T\\T\/8. 

Proof. First pick independent random local encoding 
vectors for all edges in the network simultaneously. 
Then pick an error pattern F 6 T. For this F, execute 



Steps 1 and 2 in page 3 and compute the global encod- 
ing vectors 6 F (e)'s belonging to Fg + ' F '. For each cut 
in the network, test whether Bf's are linearly indepen- 
dent for all t. This test fails with probability at most 
\T\/q by the proof of [6, Theorem 9] provided that this 
tests succeed on all the upstream cuts and n > k + 2a. 

In the proposed algorithm in Fig. 2, we test linear 
independence of Bf's on \E\ cuts in Step (**), which is 
sufficient to garuantee the decodability of the informa- 
tion by every sink. By the same reason, for each sink 
to be able to correct errors in F, one needs to consider 
linear independence only on at most \E\ such cuts with 
random choice of local encoding vectors. By the union 
bound, the probability that the the independence tests 
fails for any of |T| sinks in any of the \E\ cuts in any of 
the \T\ error patters is at most 5 if q > \E\\T\ \T\/6. □ 

Jaggi et al. [5] do not provide an estimate on the 
relation between the success probability of their algo- 
rithm and the field size q. Their method [5] uses coding 
among different time intervals and thus introduces de- 
lays while our methods do not introduce extra delay. 
In addition to this, a-error correcting codes by con- 
structed by the proposed methods allow sinks to cor- 
rect less than a errors, while the method in [5] does not. 
The advantage of the method in [5] over the proposed 
methods in this paper is that their method allows effi- 
cient decoding of information by every sink, while our 
proposed methods require exhaustive search of trans- 
mitted information. 

We summarize the comparison among the pro- 
posed algorithms and [5], [11] in Table 1. 

4.3 Relation to the robust network coding 

We clarify the difference between the robust network 
coding in [7, Sect. V],[6, Sect. VI] and the network 
error-correcting codes with known locations of errors 
[12]. A network error correcting codes that can correct 
errors on a known locations F C E is a robust network 
coding tolerating edge failures on F. However, the con- 
verse is not always true. Consider the network consists 
of three nodes {s,t, v} with two directed edges from s 
to v and one directed edge from v to t. The source is 



MATSUMOTO: OPTIMAL NETWORK ERROR-CORRECTING CODES 



7 



s and the sink is t. The intermediate node v sends to 
t the sum of two inputs from s. This network coding 
tolerate single edge failure between s and v but cannot 
correct single error between s and v. 

5. Concluding remarks 

In this paper, we proposed an algorithm constructing 
network error-correcting codes attaining the Singleton 
bound, and clarified its relation to the robust network 
coding [6, Sect. VI]. 

There are several research problems that have not 
been addressed in this paper. Firstly, the proposed de- 
terministic algorithm requires tests of linear indepen- 
dence against ('^') sets consisting of k + mf vectors, 
which is really time consuming. It is desirable to have 
a more efficient deterministic construction algorithm. 

Secondly, since there seems no structure in the con- 
structed code, the decoding of the original information 
at a sink t requires the exhaustive search by t for pos- 
sible information from the source and possible errors. 
It is desirable to have a code with structure that allows 
efficient decoding. 

Finally, the case \T\ — 1 and \E\ — n includes the 
ordinary error correcting codes as a special case. Sub- 
stituting \T\ = 1, \E\ = n and 2a = n — k into Eq. (1) 
gives q > ( n " fe ), which can be regarded as a sufficient 
condition for the existence of the MDS linear code. On 
the other hand, a well-known sufficient condition for the 
existence of the MDS linear code is q > n — 2, which 
suggests that Eq. (1) is loose and that there is a room 
for improvement in Eq. (1). 
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